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Abstract 

Diamagnetically coupled magnetic field can be used as a probe of nodal positions in unconventional superconductors. 
The heat capacity depends on the angle (jto between the magnetic field and the nodal directions. We show that the 
anisotropy C{(f>o) persists even in systems with strong paramagnetic coupling to the electrons's spins. 
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Anisotropy of thermodynamic properties of a su- 
perconductor in rotated magnetic field is often used 
to determine the location of the nodes in the su- 
perconducting gap at the Fermi surface. Magnetic 
field serves as a directional probe since the circulat- 
ing screening supercurrents around vortices excite 
quasiparticles differently, depending on the orienta- 
tion of the field with respect to the nodes. While 
initial work suggested that the density of states and 
the specific heat have a minimum when the field is 
applied along a node, we recently argued that the 
situation is more complex, and minima and maxima 
interchange in the T-H phase diagram [1,2]. 

Several superconducting materials, including 
heavy- fermion CeCoIns, have strong suppression of 
superconductivity due to paramagnetic coupling of 
electron spins to the field. This effect is insensitive 
to the field orientation and competes with the or- 
bital coupling. We investigate how the anisotropy 
of the heat capacity is affected by this competition. 

We follow the quasiclassical method [2,3] and 
include interaction of electron spins with the field 
via Zeeman term. The 4 x 4(spin-up-down-particle- 
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hole) Green's function g(R, p; 
berger equation [5] 



satisfies Eilen- 



[ers ~ Vorb - vz 



,g]+^v/• Vff = 0.(l) 



Here A(R, p) is the mean-field order parameter, 
(R; e) is the impurity self-energy. The orbital 
coupling of magnetic field is via the vector poten- 
tial A(R), while the Zeeman term couples B to the 
magnetic moments of electrons ^ ~ (g/2)fiB, 



Vorb = --V/A T3 

c 




In the absence of spin-orbit interaction we choose 
the direction of the field to be the spin quantization 
axis. Then all matrices split into independent blocks, 
for spin up and down. The equations for the two spin 
directions are independent, and differ only by a spin- 
dependent energy shift. For example, the equation 
for the off-diagonal components of ^(R, p; e) is 



-2z(eTMS)+v/(p) (^Vr-?:^A(R) 

/«^(p,R;e) = 2Aigf^^(p,R;£). (2) 
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Fig. 1. The heat capacity anisotropy for Z = 0.4 at differ- 
ent temperatures. Anisotropy curves correspond to points 
indicated in the right panel of Fig. 2. (The upper curves are 
shifted down to be on the same scale.) 

Here the energy e = £ — and order param- 
eter A = A + Airnp.i.i are renormalized by im- 
purity self-energies. In analogy to Ref.[2], we take 
the d^2_y2 gap A = A(R)V2cos20, with spatial 
structure given by Abrikosov vortex lattice, and do 
not consider a possible additional modulation due 
to Fulde-Ferrell-Larkin-Ovchinnikov state. We solve 
these equations for Green's function using Brandt- 
Pesch-Tewordt approximation. [4] To these we add 
self-consistency equations on aimp and A. The lat- 
ter is the only place where two spins add, A(R, p) = 
;^tanh^/dp'y(p,p')i(/T + /i). 

The specific heat is calculated from the density of 
states for two spins N^-^i using. 



de 



iN^{T,cj,o;e)+N^{T,(bQ-,e)) 
4T2 cosh^(e/2r) 



which is valid at low temperature, T ^ Tc- The 
calculations arc done for a quasi-cylindrical Fermi 



surface, 



Pf 



pI 



Pv 



pjcos{2spz/r Pf), with 



parameters r = s = 0.5 that ensure the 3D na- 
ture of the vortices. [2] Parameter Z ~ fiBo/2TTTc, 
where Bq = (c/i/2e)/27r^Q, and the in-plane co- 
herence length ^0 — characterises the 
strength of the Zecman term. 

Figure 1 shows the anisotropy of the specific heat 
as a function of T and B. For B in the nodal direc- 
tion, 00 = 45°, C{(f>o) has minima or maxima de- 
pending on the temperature and field range. This is 
similar to the case of purely orbital coupling, Z = 0, 
and is due to the interplay between excitation and 
scattering of quasiparticles by magnetic field, which 
is affected by the quasiparticles energy (~ T), and 
magnitude and the orientation of the field. [2] 



Fig. 2. Phase diagram of the C'{T, tpo) anisotropy for differ- 
ent strength of the paramagnetic term. In the shaded areas 
C{<j>o) has a minimum when the field is along a nodal di- 
rection. The arrows in the right panel indicate onset of the 
first order transition occurring for Z* > 0.35. 

Fig. 2 shows regions of minima and maxima at 
00 = 45° in the T-B plane, and their evolution 
with increasing Zeeman term Z . Of the two areas 
where C{(f>o) has a minimum at 45°, the low-T, low- 
B regime is only weakly affected by variation in Z. 
In contrast, the region at high T expands as we in- 
crease Z, at the expense of the area where a max- 
imum at (j>o = 45° is observed. Note that the in- 
plane anisotropy of Bc2 (present for Z = 0) quickly 
disappears as we turn on the strength of the Pauli 
coupling. 

The anisotropy of C((/)o) is usually measured be- 
low Tc/2 [6,7], where the effect of the Zeeman cou- 
pling is weak. While increased Z expands the high-T 
"minimum" region, the amplitude of the anisotropy 
there is probably still below the experimental reso- 
lution. We conclude that the Zecman coupling has 
weak effect on the observed anisotropy. 
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